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Abstract We demonstrate the existence of phase fluc- 
tuations in elongated Bose-Einstein Condensates (BECs) 
" and study the dependence of those fluctuations on the 
system parameters. A strong dependence on tempera- 
\ ture, atom number, and trapping geometry is observed. 

Phase fluctuations directly affect the coherence proper- 
[ ties of BECs. In particular, we observe instances where 
the phase coherence length is significantly smaller than 
the condensate size. Our method of detecting phase fluc- 
tuations is based on their transformation into density 
modulations after ballistic expansion. An analytic the- 
■ ory describing this transformation is developed. 
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■ 1 Introduction 

_ _ i 

Bose-Einstein condensates of weakly interacting gases, 
such as alkali atom vapours, constitute very well control- 
lable macroscopic quantum systems. For extremely low 
temperatures (T — > 0) the condensate is well described 
by an effective macroscopic single particle wavefunction, 
occupied by millions of atoms Q. The corresponding 
macroscopic phase is related to many fascinating prop- 
erties of BEC such as its coherence, superfluidity, and 
effects known from superconductivity as the Josephson 
effect. In particular, the coherence properties are essen- 
tial for very promising applications of BEC such as mat- 
ter wave interferometry [|],[|||,||,1|,[E||] or atom l asers 
H which rely on BEC as a source of coherent matter 
waves. The phase coherence of condensates was shown 
by imaging the interference pattern of two independent 
condensates which were brought to overlap (|]. Also us- 
ing an interference technique it has been found that a 
trapped BEC has a uniform spatial phase H and there- 
fore the coherence length is just limited by the size of 
the condensate. This result has also been obtained by 



measuring the momentum distribution in the radial di- 
rection of a cigar shaped condensate by a spectroscopic 
technique pQ| . These experiments have focused on the 
coherence properties of almost pure condensates, pre- 
pared at temperatures well below the critical tempera- 
ture. The temperature dependence of the coherence of 
a BEC was studied in a double slit' experiment where 
an atom laser beam was extracted from a cigar shaped 
BEC at different radial positions Q. A reduction of in- 
terference fringe visibility was observed for increasing 
temperature with the interference pattern being repro- 
ducible and thus indicating that the relative phase of the 
condensate fraction was not fluctuating randomly. 

Even though the uniform spatial phase of a conden- 
sate has been confirmed in several experiments, it is not 
an obvious or even general property of BEC at finite 
temperature. The 'fragmentation' into independent con- 
densates with randomly fluctuating relative phase was 
often discussed in the context of the nucleation process 
of a condensate jy]. There, the equilibrium state of the 
system was assumed to be a 'pure' condensate without 
phase fluctuations. However, before the condensate is 
completely formed phase fluctuations between different 
regions of the condensate are expected even though den- 
sity fluctuations are suppressed. But also for the equi- 
librium state of a quantum system it is expected that 
low-dimensional (ID and 2D) quantum gases differ qual- 
itatively from the 3D case with res pec t to statistical and 
phase correlation properties |]l2l, p^|]bl|l . Low-dimensional 
quantum degenerate gases, where the trap energy-level 
spacing is larger than the interaction energy between 
atoms, have been experimentally realized in atomic hy- 
drogen [|l5|, Sodium |HJ, and Lithium jl7j but phase 
correlation properties have not yet been experimentally 
investigated. Recently, it was shown theoretically [[f8) 
that for very elongated condensates phase fluctuations 
can be pronounced already in the equilibrium state of 
the usual 3D ensemble, where density fluctuations are 
suppressed. This was experimentally and theoretically 
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investigated in JlSf] . The phase coherence length in this 
case can be much smaller than the axial size of the sam- 
ple. 

In this paper we present detailed experimental and 
theoretical studies of phase fluctuations in elongated BEC. 
The dependence on the trapping geometry, temperature, 
and the number of condensed atoms is investigated. We 
demonstrate that the spatial phase of the condensates 
undergoes random fluctuations with an average value 
determined by the system parameters. In particular, our 
measurements show that the average phase fluctuations 
increase with tighter radial confinement. Thus, they are 
especially pronounced for very elongated geometries as 
are used to study the transition from 3D to ID degener- 
ate quantum gases ]l^ , ^9[ . 

BEC has recently been achieved in elongated micro- 
circuit geometries in which radial trapping frequencies of 
tens of KHz or even MHz are possible [g0| . Since this is 
more than an order of magnitude bigger than the radial 
frequencies at which we have observed phase fluctua- 
tions our results should be especially relevant for these 
systems. Phase fluctuations should also be considered if 
tight confinement waveguides are used for BECs || or 
guided atom lasers beams, even though our results di- 
rectly apply only to the equilibrium state of the system. 

Note, that phase fluctuations lead to a broadening of 
the momentum distribution, which is the Fourier trans- 
form of the single particle correlation function, defining 
the coherence length. Therefore, only in the absence of 
phase fluctuations, a BEC constitutes a state where ulti- 
mate control over the motion and position of atoms, lim- 
ited only by Heisenberg's uncertainty relation, is achieved. 

Fluctuations of the phase of a Bose condensate are 
due to thermal excitations and always appear at finite 
temperature. Their experimental characterization con- 
stitutes a test of many-body theories at finite temper- 
ature. If the wavelength of these excitations is smaller 
than all dimensions of the condensate they have a 3D 
character and do not lead to pronounced phase fluctua- 
tions. In elongated condensates however, low energy ax- 
ial excitations have wavelengths larger than the radial 
size of the sample and therefore acquire a pronounced ID 
behaviour leading to axial phase fluctuations, although 
density fluctuations of the equilibrium state are still sup- 
pressed. The coherence properties can be significantly 
altered as compared with previous observations. In par- 
ticular, the coherence length, i.e. the distance at which 
the single particle correlation function falls to e , can 
be much smaller than the axial size of the condensate. 
This is not in contradiction to previous coherence mea- 
surements since they were performed in rather spherical 
traps H , in the radial direction of cigar shaped conden- 
sates and at low temperature p|,^0|. 

Our method to study phase fluctuations is based on 
ballistic expansion. We show that the original phase dis- 
tribution is mapped into the density distribution during 
time-of-flight. The density modulations after ballistic ex- 



pansion are a direct measure of the original phase fluc- 
tuations of the trapped condensate. 

The paper is organized as follows. It starts with a 
theoretical discussion of phase fluctuations including a 
detailed explanation of our analytic theory |n| for the 
appearance of density modulations in the ballistic expan- 
sion. Then we present our experimental studies demon- 
strating the existence of phase fluctuations in BECs and 
investigating their dependence on system parameters. 
Finally, alternative methods to study phase fluctuations 
are discussed, including simulations of the interference 
of phase fluctuating condensates. 



2 Theoretical description 

2.1 Phase fluctuating condensates in elongated 3D traps 

In a standard equilibrium situation in 3D traps the fluc- 
tuations of density and phase of the Bose-Einstein con- 
densate are only important in a narrow temperature 
range near the BEC transition temperature T c (criti- 
cal fluctuations). Outside this region, the fluctuations 
are suppressed and the condensate is phase coherent. 
This picture precludes the interesting physics of phase- 
fluctuating condensates, which is present in 2D and ID 
systems (see 12>13f and refs. therein). 

In the Thomas- Fermi regime (see e.g. [Q), when the 



nonlinear mean field interaction energy dominates the 
kinetic energy, the situation can change for elongated 
BECs . While density fluctuations in equilibrium will 
remain suppressed due to their energetic cost, the situ- 
ation can be different for the phase fluctuations. In par- 
ticular, the axial phase fluctuations can manifest them- 
selves even at temperatures far below T c . Then, as the 
density fluctuations are suppressed, the equilibrium state 
of the system becomes a condensate with fluctuating phase 
(quasicondensate) similar to that in ID trapped gases 
|12| . Decreasing T gradually reduces the phase fluctua- 
tions. 



2.2 Description of the fluctuating phase 

In this section we describe the phase fluctuations along 
the lines of Ref. Q . Let us consider a BEC at T = 
in the Thomas-Fermi regime, where the mean-field (re- 
pulsive) interparticle interaction greatly exceeds the ra- 
dial (uj p ) and axial (u) x ) trap frequencies. The density 
profile of the zero-temperature condensate has the well- 
known shape no(p, x) — n$ m (l — p 2 j ' R 2 — x 2 / L 2 ), where 
nom = p/g is the maximum condensate density, with 
p, being the chemical potential, g — 4:Trh 2 a/m, m the 
atom mass, and a > the scattering length. Under the 
condition lo p 3> lo x , the radial size of the condensate, 
R = (2/i/mWp) 1 / 2 , is much smaller than the axial size 
L = (2 M /m^)l/2. 
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The phase fluctuations can be described by solving 
the Bogoliubov-de Gennes equations |nj describing ele- 
mentary excitations of the condensate. One can write the 
total field operator of atoms as ?/>(r) = ^Jn Q (Y) exp(ic/>(r)), 
where 4>(r) is the operator of the phase, and the density 
fluctuations have been already neglected following the 
arguments discussed above. The single-particle correla- 
tion function is then expressed through the mean square 
fluctuations of the phase (see, e.g. p2|): 



(^(iW)} = Vn (r)no(r')ex P {-([^(r ) r')] 2 )/2}, (1) 



with 8<f>(r,r') — </>(r) — 
by (see, e.g., §f|) 

(j>{r) = [4n (r) 



(r'). The operator 4>(r) is given 



h.c. 



(2) 



where a,j is the annihilation operator of the excitation 
with quantum number(s) j and energy tj, fj = Uj + vj, 
and the it, v functions of the excitations are determined 
by the Bogoliubov-de Gennes equations. 

The "low-energy" axial excitations (with energies €j < 
huip) have wavelengths larger than R and exhibit a pro- 
nounced ID behavior. Hence, one expects that these ex- 
citations give the most important contribution to the 
long-wave axial fluctuations of the phase. The solution of 
the Bogoliubov-de Gennes equations for such low-energy 
axial modes gives the spectrum ej — fruj x y/j (j + 3) /4 
|plf , where j is a positive integer. The wavefunctions fj~ 
of these modes have the form 



f + M - / (j + 2)(2j + 3)gn (r) (1>1) ,x\ 



where pi 1 ' 1 ' are Jacobi polynomials. Note that the con- 
tribution of the low-energy axial excitations to the phase 
operator (H) is independent of the radial coordinate p. 



2.3 Numerical simulations 

In order to simulate numerically the effect of phase fluc- 
tuations during the ballistic expansion of the conden- 



sate we replace the operators a, and a\ in Eq. 



by 



complex Gaussian random variables ctj and a* , with the 
correlation (aj-a*,) = SjjiNj, where Nj is the occupation 
number for the quasiparticle mode j for a given (small) 
temperature T. Such a random phase reproduces cor- 
rectly the phase correlations, which for not too large 
\x — x'\/L < 0.4, behave as 



([64>(x,x')f) T = Sl\x 
where the quantity 8\ is given by 



AIL, 



6 2 L (T)=32nT/15N {hu x ) 2 , 



(4) 



(5) 



and I ' a2 
length. It is the distance at which the phase factor of the 



L/8\ can be interpreted as a phase coherence 



2tt 
In 
Ojc 

-l7t 

-2jc- 



In 
In 

,_, Ore 
1 "If 



2tz 
In 
On 
-In 
-In 




-0.4 



0.0 

x/L 



0.4 0.8 



Fig. 1 Typical phase patterns for three different trap as- 
pect ratios: A = 10 (top), A = 100 (middle), and A = 1000 
(bottom). For all cases lu x = 2n x 14 Hz, T = 0.6 T c , and 
N = 2 x 10 5 . 



single-particle correlation function (Eq. [I]) falls to its 1/e 
value. The condition l^/L = 1 determines then a char- 
acteristic temperature — 15(Tilu x ) 2 No /32/j,, where iVo 
is the number of condensed atoms. For < T c , which 
is the case for most of our measurements, one expects 
the regime of quasicondensation for the initial cloud in 
the temperature interval < T < T c Q. 

Figure [l] shows examples of simulated phase distribu- 
tions for various trap geometries parameterized by the 
ratio of trapping frequencies A = u) p /uj x . Whereas for a 
trap with A = 10 the fluctuations are rather suppressed, 
the phase of the BEC fluctuates for a trap with A = 1000 
by more than 2ir. This dramatically changes the proper- 
ties of the condensate, especially in phase sensitive ex- 
periments. 

The appearance of stripes in the process of ballis- 
tic expansion can be understood qualitatively as follows. 
As mentioned above, within the equilibrium state of a 
BEC in a magnetic trap the density distribution remains 
largely unaffected even if the phase fluctuates Jn| . The 
reason is that the mean-field interparticle interaction 
prevents the transformation of local velocity fields pro- 
vided by the phase fluctuations into modulations of the 
density. However, after switching off the external trap, 
the mean-field interaction rapidly decreases and the axial 
velocity field is then converted into a particular density 
distribution. We have performed numerical simulations 
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of the 3D Gross-Pitaevskii equation (GPE) to under- 
stand quantitatively how phase fluctuations lead to the 
build up of stripes in the density distribution. We as- 
sumed that initially (just before opening of the trap) the 
condensate had an equilibrium density profile, which has 
been calculated using the standard imaginary evolution 
of the GP equation. We have imprinted on it a random 
fluctuating phase <fr(x), as described above, and evolved 
the condensate in free space. 



2.4 Analytic results 

Alternatively, the appearance of the stripes and their 
statistical properties can be described analytically using 
the self-similar solutions of the GP equation valid for the 
expanding cloud in the Thomas-Fermi regime (see e.g. 
pi| ) . In the following we shall assume for simplicity that 
the condensate is an infinite cylinder. This assumption is 
justified since the typical size of the excitations is much 
smaller than the axial size of the condensate. Therefore, 
at the end of the calculation, the unperturbed density 
will be substituted by the corresponding local density. 
In this way, the condensate without initial fluctuations 
evolve according to the self-similar solution 



(6) 



where b 2 p (t) = 1 + Lu 2 t 2 , (j) — ^--jf-p 2 , p is the radial 
coordinate, and n is the Thomas-Fermi density profile. 
The equations which determine the ballistic expansion in 
presence of fluctuations can be obtained by linearizing 
around the self-similar solution, for the density no + Sn, 
and the phase 0o + 0. Introducing this expressions in the 
corresponding Gross-Pitaevskii equation, we obtain: 



21 ), we end up with the equation: 



Sn 



2mb 2 p 



VtSn- 



4m 2 



-VtSn = 0. 



(10) 



Expanding into the Bogoliubov eigenmodes of the sys- 
tem, then 

e 2 k {t)Sn k = (11) 



h 2 Sn k 



where 



6fc(t) 



\ ph 2 k 2 i fi 4 fc 4 

2mb 2 



4m- 



(12) 



is the Bogoliubov spectrum. 

For short times uj p t < 1 the phonon part of the spec- 
trum is dominant, i.e. e k 

(*) - The resultin s 

equation can be solved in terms of hypergeometric func- 
tions. For larger times, the free particle part of the spec- 
trum dominates, i.e. e k = 2 , and the equation can be 
solved in terms of Bessel functions. The two regimes may 
be matched, using the asymptotic expansions of both hy- 
pergeometric and Bessel functions. In this way, we end 
up with the analytic expression for the relative density 
fluctuations 



Sn 



— =2rr( e '« sin 
n ^ 



phiOp 



(13) 



where the sum extends over the axial modes ej — Tilo x y/ j(j + 3)/4, 
r = ujpt, and <j)j is the contribution of the j— th mode to 
the phase operator in Eq. (||). 

By substituting the operators by the corresponding 
c-numbers as discussed above, we have checked that the 
analytical expressions agree very well with the numerical 
results (see also Fig. ||). 



-^Ocf) , 

6p m 



n (, 



gngSn 



± 6 (Sn 
46* 



no 



hV 2 Jr 
4m 



(7) 
(8) 



where O = -{%/m)(V p m - V 
By combining Eqs. (0) and 

gnoV 2 8n h 2 V^Sn 



n V 2 ,), and p' = p/b p 



we obtain: 



mb 2 



Am 2 




Ab 2 m x no 



(9) 

The last term at the rhs of Eq. (S) can be neglected 
if [ijhiOp 3> 1 (i.e. in the Thomas-Fermi limit). The 
third term in the lhs of the equation (quantum pres- 
sure term) is for short times smaller than the second 
one on the lhs, but becomes comparable with it for times 

t ~ (\/uj p )\J /im/h 2 k 2 . This becomes an important point 
as discussed below. Averaging over the radial profile (em- 
ploying the expansion in powers of p as discussed in Ref. 



2.5 Averages 

From Eq.Q one obtains a closed relation for the mean 
square density fluctuations a 2 by averaging (Sn/no) 2 
over different realizations of the initial phase. Taking 
into account that {(a + a*)(a > + a|,))/4 = Nj8jj</2, 
then the mean square fluctuations are given by: 



( 5n(x,t) y\ 
\\n (x,t)J j 

where 

C(N ,uJp,uJ x ,x,t) 2 



T 



C(A^o,^ P ,^,x,t) 2 , (14) 



1 ^ 

iVsin 2 

2 ^ 



(J+3/2) 2 
4q(1- {x/LY) 



■ (j + 3/2) 2 



ln(2u p t) 



(j + 2)(2j + 3) x (1 ,d 
with a = p,/huj 2 t. 



©> 2 



(15) 
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Fig. 2 Typical density profile after 25 ms time-of-flight. 
Top: numerical simulation without taking into account a 
limited experimental resolution. Bottom: numerical simula- 
tion (dotted line) compared with the analytic theory (solid 
line), both taking into account an experimental resolution of 
a = 3/im. All profiles were calculated for uj x — 2ir x 14 Hz, 
Up = 2tt x 508 Hz, No = 2 x 10 5 , and T = 0.5 T c . For all the 
figures the same initial phase pattern was used. 

Using the quasiclassical average (P^'^^x)) 2 ~ 4(1 — 
x 2 )~ 3 / 2 /nj and transforming the sum over j into an in- 
tegral, for the central part of the cloud (x « 0) we obtain 




Note that Eq. ([l6]) provides a direct relation between 
the observed density fluctuations and temperature, and 
thus can be used for thermometry at very low T. 

This closed expression would provide an accurate de- 
scription of the density fluctuations in the absence of any 
experimental limitation. However, in practice, the obser- 
vation of the density fluctuations is limited by the spa- 
tial resolution of the experiment. This fact can be easily 
taken into account in the calculations by substituting 
Pj 1,1 '(x/L) in Eq. JIB] ) by the corresponding function 



convoluted by a resolution function £(x) = e~ x l G /y/na, 
where a characterizes the experimental resolution (Fig. |^). 

Other experimental limitations can be easily taken 
into account by following a similar procedure. In partic- 
ular, a smoothing of the observed density ripples is pro- 
duced if the laser that integrates the condensate column 
density in the absorption detection scheme is not exactly 
parallel to the density stripes. This effect can be easily 
incorporated into the calculations, by filtering the den- 
sity distribution with a cut-off function in the Fourier 
spectrum of the form exp(ikRb p (t)6) sm(kRb p (t)6)/kTb p (t)8, 
where 6 is the angle between the laser and the density 
ripples. We return to this point in the experimental sec- 
tion. 



3 Experimental Results 

3.1 Experimental Setup 

The experiment was performed with Bose-Einstein con- 
densates of 87 Rb atoms in the \F = 2, nip = +2) hy- 
perfine ground state. As described previously |2R^6| we 
load a magneto-optical trap with a few times lCr atoms 
from a chirp slowed thermal beam. This is followed by 
a short period of subdoppler cooling and optical pump- 
ing into the desired magnetic sublevel. The atomic cloud 
is then loaded into an Ioffc-Pritchard type (cloverleaf) 
magnetic trap and finally adiabatically compressed to 
allow efficient rf evaporative cooling. The fundamental 
frequencies of our magnetic trap are cu x = 2ir x 14 Hz 
and u>p = 2tt x 365 Hz along the axial and radial direc- 
tion, respectively. Due to the highly anisotropic confin- 
ing potential with an aspect ratio A = lo p /lo x of 26 the 
condensates are already elongated along the horizontal 
x-axis. 

In order to study the dependence of phase fluctua- 
tions on the trapping geometry we realized a wide range 
of radial confinement strengths. To allow for stronger 
radial confinement we used a holographically generated 
blue detuned Laguerre-Gauss mode (TEM^) laser beam 
to form a combined magnetic and optical dipole poten- 
tial trap H . In this combined trap BECs were produced 
in a two step evaporation procedure. The atoms were 
first cooled in the 'pure' magnetic trap to a temperature 
slightly above the transition temperature T c . The optical 
dipole potential was then turned on adiabatically, and fi- 
nally the desired temperature was achieved by rf evapo- 
ration in the combined potential. On the other hand less 
elongated condensates were produced by using a higher 
offset-field for the magnetic trap leading to a weaker ra- 
dial confinement. 

Our measurements were performed for an axial trap 
frequency of u> x = 2ir x 14 Hz and radial frequencies uj p 
between 27r x 138 Hz and 27r x 715 Hz corresponding to 
aspect ratios A between 10 and 51. After rf evaporative 
cooling to the desired temperature, we wait for 1 sec 
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Fig. 3 (a)-(c): Absorption images for various times-of-fiight 
and temperatures with u) x = 2nx 14 Hz and u p — 2n x 365 Hz. 
(d): Density profiles for the clouds displayed in (c) integrated 
along the radial direction. In the case of T < 200 nK no ther- 
mal component was visible and the temperature was esti- 
mated to be < T c /2. 



(with rf 'shielding') to allow the system to reach an equi- 
librium state. We then switch off the trapping potential 
within 200 fis and wait for a variable time-of-flight be- 
fore detecting the atomic cloud by resonant absorption 
imaging with the imaging axis perpendicular to the long 
axis of the trap. 



3.2 Ballistic expansion measurements 

Figure || shows typical images of the ballistically ex- 
panded clouds for various times-of-flight t and temper- 
atures T < T c . The usual anisotropic expansion of the 
condensate due to the elongated trapping geometry is 
clearly visible. The line density profiles reflect the parabolic 
shape of the BEC density distribution. As predicted by 
the theory, we also observe pronounced stripes in the 
density distribution. On average these stripes are more 
pronounced for high temperatures (right column of Fig. ||) 
and long times-of-flight [Fig. ||(c)] indicating the build- 
up of the stripes during the ballistic expansion. 



3.3 Experimental characterization of phase fluctuations 

In order to determine the dependence of phase fluctua- 
tions on the experimental conditions we study the for- 
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Fig. 4 Distribution of phase fluctuating condensates as a 
function of temperature and condensate number for four dif- 
ferent trap geometries. White points indicate the absence of 
detectable structures (<t B ec < 1-5ct t ), black points the exis- 
tence of significant structures (<7 B ec > 2<r T ). 



mation of the density modulations as a function of tem- 
perature and trapping geometry. All measurements were 
performed for times-of-flight between 18 ms and 25 ms. 
For a quantitative analysis we compare the observed den- 
sity distribution with the Thomas-Fermi distribution ex- 
pected for a condensate without fluctuations. For each 
image the observed density distribution was integrated 
along the radial direction [Fig. |^(d)] and then fitted by a 
bimodal function with the integrated parabolic Thomas- 
Fermi density distribution and a Gaussian for the ther- 
mal cloud. In Fig ||(d) we show the integrated density 
distribution close to the BEC region, thus the thermal 
component is hardly visible. The standard deviations of 
the experimental data from the fit were calculated in the 
central region of the BEC (half width of full size) , <r BEC , 
and in the thermal wings, a T . The standard deviation in 
the thermal wings characterizes our detection noise. To 
account for particle number changes the standard devi- 
ations were normalized to the fitted peak density uq. 

Figure |4| shows the distribution of phase fluctuating 
condensates as a function of temperature and condensate 
number for various trap geometries. In this figure white 
points indicate the absence of detectable structures in 
the density distribution (<r BEC < 1.5<t t ), whereas black 
points indicate that significant structures larger than the 
experimental noise level (er BEC > 2cr T ) were observed, 
i.e. the existence of phase fluctuations could be clearly 
detected. The temperature and particle number of each 
condensate were determined by 2D fits to the absorption 
images. The temperature was determined from the width 
of a Gaussian distribution fitted to the thermal wings, 
the corresponding condensate number from the integral 
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Fig. 5 Measurement of ctbec /no versus temperature in four 
different trap geometries. The dotted lines represent the av- 
erage detection noise ar/no- All data was taken for T < T c . 
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Fig. 6 Average standard deviation of the measured line 
densities [ (obec/^o) 2 ]tlp as a function of temperature. The 
temperature-axis is scaled according to the expected depen- 
dence on the other experimental parameters (Eq. |l4|). The 
dashed line is a fit of a square root function to the data re- 
flecting the theoretical expected dependence. 



over the Thomas- Fermi part of a bimodal fit. The sta- 
tistical uncertainty in the temperature determination is 
typically 15 % and less than 10 % for the particle num- 
bers. 

It is clearly visible in Fig. |] that the number of 
condensates showing detectable phase fluctuations is in- 
creasing with the aspect ratio. Furthermore, the num- 
ber of realizations without detectable phase fluctuations 
is growing for low temperatures and high particle num- 
bers. Finally, there is no clear transition line between the 
two regimes. There is rather a broad region in which we 
observe both, condensates with and without detectable 
phase fluctuations indicating the statistical character of 
the phase fluctuations. Whereas for the weakest radial 
confinement with A = 10, no significant structures were 
observed outside the close vicinity of the critical temper- 
ature, we detected pronounced phase fluctuations in a 
broad temperature range in the case of our tightest trap 
A = 51 [Fig. |4j] . According to Eq. ||phase fluctuations can 
be reduced below any detection level for sufficiently low 
temperatures. However, for high aspect ratio traps this 
requires very low temperatures at high particle numbers 
in the condensate, making it experimentally difficult to 
access. For all traps the BECs were produced by evap- 
orating in the final potential except in the case of the 
A = 10 trap. This trap was realized by evaporating in the 
tighter A = 26 trap and then adiabatically reducing the 
radial confinement. Thus the measurements show that 
by changing the trapping potential adiabatically we are 
able to decrease the amount of phase fluctuations. The 
dynamics of the appearance and disappearance of the 
phase fluctuations remains to be studied systematically 
in future work. 



Fi gure |^ shows the temperature dependence of (7 BEC 
for various trap geometries. Since the initial phase of a 
Bose condensate is mapped into its density distribution 
after ballistic expansion, the quantity <7 BE c is a direct 
measure of the initial phase fluctuations. Note however, 
that this method reflects the instantaneous phase of the 
BEC at the time of release and therefore images taken at 
the same initial conditions can look significantly differ- 
ent. Indeed, we observe a large spread of our experimen- 
tal data (see Fig. ||). The scatter in these cases is mainly 
due to the statistical character of the phase fluctuations 
as well as due to the uncertainties in the temperature de- 
termination. For all traps the highest temperature data 
correspond to values close to T c (T as 0.9T C ) which is in- 
creasing for tighter confinement. On average the phase 
fluctuations continuously decrease with falling temper- 
ature and get more pronounced with increasing aspect 
ratio. The reduction of the observed phase fluctuations 
for lower temperatures is due to both, the reduced exci- 
tation spectrum at lower temperature and the increas- 
ing number of condensed atoms. Due to the loss of par- 
ticles caused by evaporation, lowering the temperature 
reduces the total number of particles but the fraction of 
condensed atoms relative to the total particle number is 
increased as shown in Fig. [|. It is not possible to de- 
termine a cut-off for the phase fluctuations, rather they 
decrease until they cannot be resolved below our noise 
limit. Hence all experiments with BECs at finite tem- 
perature in tightly confining elongated potentials will 
be subject to axial phase fluctuations. 
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3.4 Evaluation of averaged phase fluctuations 

Due to the statistical character of phase fluctuations ev- 
ery experimental realization has a different phase distri- 
bution. Therefore, we average the observed standard de- 
viations for a small range of experimental parameters to 
obtain general information about the phase fluctuations 
and to compare our measurements to the theoretical pre- 
diction. According to Eq. ( [Til) the standard deviation 
of the fluctuations is expected to depend on the square 
root of the temperature. This behaviour can not be ob- 
served in Fig. |B| since the condensate number changes 
with temperature. Therefore, we use that equation to 
scale out the influence of condensate number, trapping 
potential, and time-of-flight. In figure ^| the measured 
standard deviations are plotted as function of T/T^C 2 , 
where T$ and C (C taken at x=0) are calculated from 
the measured number of condensed atoms, trapping fre- 
quencies, and expansion times. The figure shows that 
all data points fall on an universal square root shaped 
curve. This means, that our experimental results follow 
the expected dependence very well. 

As a main result, with the direct link of the phase 
fluctuations in the magnetic trap to the observed den- 
sity modulations given by Eq. (|l6|), our measurements 
confirm the predicted behaviour of phase fluctuations in 
elongated BECs. 

As explained above, a precise quantitative compari- 
son needs to take into account a reduced contrast due to 
the limited experimental resolution and a possible tilt of 
the detection laser beam with respect to the structures. 
If we account for a resolution of our imaging system of 
a = 3/im, the predicted standard deviations are about a 
factor of 2 bigger than the observed deviations. A pos- 
sible reason might be a small tilt of the detection laser 
beam. For our parameters a tilt of only 4° reduces the 
observed modulations by a factor of approximately 2. 

Note, that our measurements were performed in elon- 
gated geometries but our condensates still were in a 
3D-regime in the sense that the chemical potential was 
greater than the transverse level spacing. Nevertheless, 
most of our measurements, which exhibit fluctuations 
well above our noise level, correspond to the regime of 
quasicondensation in which the phase coherence length 
is smaller than the condensate size. For instance, for 
A = 51, T = 0.5 T c , and N a = 3 x 10 4 , one obtains 
ji = 3ATiuj p and 1$ ~ L/3. 



4 Discussion of alternative methods 

In principle, phase fluctuations may be observed by all 
experiments on BEC which rely on its phase properties. 
Hence, various methods may be used to study them. We 
have presented a method which is based on the trans- 
fer of phase fluctuations into density modulations after 
ballistic expansion. A great advantage of our method 



is its applicability to very different trap geometries and 
that the sensitivity depends on the chosen time-of-flight 
t (see Eq. (|l6|)). In principle for longer t the method be- 
comes more and more sensitive for initial phase fluctua- 
tions. On the other hand the decreasing density reduces 
the signal-to-noise ratio leading to an optimum time-of- 
flight for given parameters. In our case a time-of-flight 
of 25 ms allows us to measure phase fluctuations of only 
w 7r/7. Moreover, time-of-flight methods constitute a 
standard experimental tool to study properties of BECs. 
Therefore, understanding the formation of stripes in the 
density distribution is of great importance. It is instruc- 
tive to discuss two of the most successful methods used 
to study the coherence properties of BEC, namely Bragg 
spectroscopy [[[o| and interferometry || , with respect to 
the measurement of phase fluctuations. 



4-1 Bragg spectroscopy 

The velocity field of a condensate is proportional to 
the gradient of its phase. Therefore, phase fluctuations 
lead to a broadening of the momentum distribution of a 
trapped condensate which can be measured using Bragg 
spectroscopy. Since in elongated condensates phase fluc- 
tuations are predominantly provided by axial excitations 
the momentum distribution along the axial direction has 
to be measured. The average momentum distribution V 
is given by the Fourier transform (FT) of the integrated 
single-particle correlation function G(r a ) = J d 3 r{^ (r)ip(r - 
27 . Using Eq. (|l|) and the definition of the phase coher- 
ence length we get 



V(p) = TT 



exp 



exp 



d 3 ry/n (r) n (r + r a ) 



I'al 

' l 4> 



(a/ n (r) n (r + r a )%l7) 



which is a convolution of the Fourier transform of the 
density distribution and that of a Gaussian with the 
width of the phase coherence length 1$. Thus, the first 
term is inversely proportional to 1$, the second inversely 
proportional to the condensate size. In the absence of 
phase fluctuations the momentum distribution is lim- 
ited by the size of the condensate, i. e. , Heisenberg- 
uncertainty limited. In the regime of quasicondensates 
where the phase coherence length is smaller than the 
condensate size, the momentum distribution is domi- 
nated by the phase coherence length, thus increasing 
the momentum width significantly with respect to the 
Heisenberg limit. However, for a trapped condensate the 
width of the Bragg resonance is not only determined by 
this momentum width (Doppler width) but also by the 
inhomogeneous density distribution, leading to a spa- 
tially inhomogeneous mean field-shift of the resonance 
and therefore broadening it. 

We will now briefly estimate the relative contribution 
of the momentum width and the mean-field width to the 
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total linewidth of the Bragg resonance. For typical pa- 
rameters of N — 10 5 , w p = 2ttx 365 Hz, w x — 2tt x 14 Hz 
the momentum width resulting from the size of the con- 
densate is <4^si ze — v^ 2 l/8 2 ^ L — 57 Hz, where k = 
2tt/X is the absolute value of the wavevector of the Bragg 
beams. For these parameters the mean field broadening 
is Av m{ = ^8/147^ = 536 Hz. Hence the total width 

of the resonance is Av = ^^ omentum + 4< f = 
539 Hz. Here, we have taken the expressions for the widths 
from Ref. fl0|| and assumed Bragg spectroscopy with 
counterpropagating beams since this leads to the high- 
est momentum resolution. For a temperature T = 0.5T C 
the phase coherence length is, according to Eq. (||), 1$ m 
1.4L w 7r/5). Though this increases the momentum 
width roughly by a factor of 2, it leads only to a to- 
tal linewidth of approximately 547 Hz. Thus, the contri- 
bution of the momentum width to the total linewidth 
in this case is less than 2%, making it very difficult 
to detect. In contrast, our measurements show that the 
method of ballistic expansion is capable of detecting even 
such small phase fluctuations. 



4-2 Interferometry 

Interferometry has been successfully used to study and 
to characterize the spatial phase of Bose condensates. 
Therefore, it is natural to consider possibilities to ex- 
tend such methods to determine the amount of phase 
fluctuations in trapped atomic gases. In analogy to clas- 
sical optics, there are several methods to use interferom- 
etry to determine the coherence properties of a conden- 
sate wave function. These include, e.g., the double slit 
experiment with outcoupled atom laser beams Band 
interference between two overlapping condensates [0] . In 
the latter case, two spatially separated condensates were 
prepared independently and their interference was mea- 
sured after a ballistic expansion. The results indicated a 
flat phase distribution of the initial wavefunction for the 
experimental parameters used in those measurements. 
In general, interferometric measurements rely on the su- 
perposition of two wavefunctions and the interpretation 
of the resulting fringe pattern. The existence of a ran- 
domly fluctuating phase as shown in Fig^ complicates 
the interpretation. 

To visualize possible experimental results, let us first 
consider the following situation. A given phase fluctu- 
ating condensate is coherently split in two equal parts 
giving one of them a small additional velocity. Absorp- 
tion images taken of the overlapping condensates will 
show interference fringes. For condensates with a con- 
stant phase the images would show regularly spaced lines. 
However due to the phase fluctuations these lines will be 
shifted according to the relative fluctuating phase, lead- 
ing to a random, spatially varying fringe spacing. Simi- 
lar to the ballistic expansion measurements these images 




Fig. 7 (a, b) Examples of independently calculated density 
distributions of phase fluctuating condensates (solid line) for 
Lu p = 2ttx715Hz, u a = 2ttx14Hz, TV = 5xl0 4 andT = 0.6T C 
compared to that of a condensate at T = (dashed line) 
for a time-of-flight of 15 ms. (c) Interference pattern of the 
superposition of the condensates shown in a) and b) with a 
relative displacement of Ax = 40^m. 

will be different from shot to shot and allow only for a 
statistical evaluation. 

So far most interferometric measurements with BECs 
have been performed after a ballistic expansion. An ad- 
ditional complication arises in these experiments, since 
both density and phase of the sample evolve during the 
free expansion time. Let us assume that two indepen- 
dent quasicondensates are brought to partial overlap af- 
ter a time-of-flight with zero relative velocity (using e. g. 
Bragg pulses || f7[H] ) ■ At this point modulations have 
appeared in the density distribution and the phase has 
evolved from its initial random pattern. Both parts now 
contribute to the formation of stripes in the interference 
pattern, making a direct interpretation difficult. 

Figure ^ shows an example of an expected absorp- 
tion image for our experimental conditions. Both wave- 
functions <?i and ^2 of the phase fluctuating conden- 
sates were calculated after a time-of-flight of 15 ms us- 
ing the GP equation. The interference pattern is given 
by |^i (x) +$2(2: + <4x)| 2 , where Ax is the displacement 
of the two quasicondensates. Fig. |?] a) and b) show the 
density distribution of the two quasicondensates (solid 
line) displaced by Ax = 40/Ltm « 21^ in comparison with 
the Thomas- Fermi distribution of a condensate at T = 
(dashed line). Note, that the structures shown there are 
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the density modulations due to the phase fluctuations 
only. The resulting interference pattern of these two in- 
dependent condensates is then displayed in Fig. c). At 
T = the interference pattern is due to the small axial 
expansion velocity, which leads to a quadratic phase. In 
the case of the quasicondensates the profile looks con- 
siderably different, i. e. the phase fluctuations lead to 
a significant change of the experimental outcome. Due 
to the difficult interpretation of the density distribution, 
interferometric methods together with a long expansion 
time are not suited very well for analysis of phase fluctua- 
tions. However, interferometric methods that work with 
a short expansion time or even with trapped conden- 
sates might be very useful for visualization of the spa- 
tial profile of phase fluctuations. Such methods might 
include, e.g. preparation of two independent quasicon- 
densates with a small vertical separation. After a short 
time-of-flight the radial expansion of the clouds leads 
to their partial superposition, and consequently interfer- 
ence in the overlap region. In the case of constant relative 
phase in the axial direction horizontal stripes appear in 
the overlap region. Due to the phase fluctuations the 
relative phase can, however, change as a function of ax- 
ial position, which is seen as bending of the interference 
fringes in the vertical direction. Therefore, the relative 
phase profile of the clouds can be imaged by measuring 
the shape of the interference fringes. 



5 Conclusion 

We have presented detailed experimental and theoreti- 
cal studies of phase fluctuations in the equilibrium state 
of BECs. A strong dependence on the trap geometry 
and temperature has been found in agreement with the 
theoretical prediction. We have shown that phase fluc- 
tuations are a general property of elongated conden- 
sates. By measuring the phase fluctuations and compar- 
ing the temperature with we have demonstrated in- 
stances, where the phase coherence length was smaller 
than the axial size of the condensate, i.e. the initial cloud 
was a quasicondensate. Our results set severe limita- 
tions on applications of BECs in interferometric mea- 
surements, and for guided atom laser beams. Our exper- 
imental method combined with the theoretical analysis 
provides a method of BEC thermometry. Further studies 
of the effect of phase fluctuations, e. g. on the superfluid 
properties, will give additional insight to the behaviour 
of degenerate quantum gases at finite temperature. 
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